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Abstract 

A leading-power factorization formula for weight functions relating the B — > X s ^f photon 
spectrum to arbitrary partial decay rates in B —* X n Z~P is derived. These weight 
functions are independent of the hadronic shape function and allow for the determination 
pJIj ■ of \V u b\ in a model-independent way. We calculate the weight function in renormalization- 

group improved perturbation theory to complete next-to-next-to leading order at the jet 
scale jti| ~ m^AQCD and to next-to leading order at the hard scale fih ~ mb- First-order 
power corrections are also included, where a model-dependence enters via the appearance 
of subleading hadronic shape functions. 
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1 Introduction 

The determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element \V u b\ from in- 
clusive semileptonic B — > X u l~v decays requires theoretical predictions for partial decay rates, 
which are then compared to their experimentally measured values. Because of a dominating 
b — > c background in a large portion of phase space, this procedure is adopted for a variety of 
restricted regions in phase space as obtained by, e.g., accepting only events with one or more 
of the following features: charged- lepton energy E\ > Eq, hadronic invariant mass Mx < Mq, 
leptonic invariant mass q 2 > qfi, hadronic P + < A. Here, P± = Ex =F \Px\, where Ex denotes 
the energy and Px the three-momentum of the final hadronic state in the 5-meson rest frame. 
For many of these cuts a hierarchy of energy scales exists in the decay process, for example 
P + ~ Aqcd <C Mx ~ a/^Aqcd "C mt, and shape-function effects become important [TJEllHj- 
The theoretical expressions for differential decay rates in this region of phase space factorize 
into hard functions at the scale fit ~ m&, and the convolution of jet functions and shape func- 
tions at the scale /ij ~ a/^Aqcd jHIEj- While the jet functions are perturbatively calculable, 
shape functions are non-perturbative objects that capture all strong-interaction effects below 
the scale /ij. At leading power the jet function J and shape function S are universal, and en- 
ter the QCD-factorization theorems for both the triple differential decay rate in B —>■ X u l~v 
decays and in the B — > X s ^ normalized photon spectrum jHJ E] . 

One strategy for the inclusive determination of \V u b\ is to use the B — > X s j photon spec- 
trum to extract the leading shape function S, which then allows for the calculation of arbitrary 
semileptonic partial decay rates jSj. It was shown in that reference that the factorization ap- 
proach can be applied to all commonly used kinematic cuts. In practice this program is realized 
by adopting a parameterizable model for the shape function, and fitting the parameters using 
information of the measured photon spectrum. It was emphasized in jS] that such a model is 
only acceptable if the values for parameters are stable when fitted to different aspects of the 
photon spectrum, such as moments of it, or its functional form (provided that the resolution 
is coarse enough to smear out hadronic resonances). With improving data on the photon 
spectrum such an approach might require further and further refinement of the models used. 
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Figure 1: Relevant part of the phase space of semileptonic B decays. The light-gray region 
denotes the portion where a background from B — > X c l~v decays exists, while the dark- 
gray area is free of this background. Events with final-state hadronic invariant mass Md are 
located on the line between these two regimes. Events with P + < M^/Mb are located below 
the dashed line. As a slightly more complicated example we also plot the borderline in phase 
space for events with Mx < 1.7 GeV and q 2 > 8 GeV 2 as a dotted line. 



A different strategy is to eliminate the necessity for the extraction of the leading shape 
function, and to use the experimental data directly. Such ideas have been investigated previ- 
ously in [3 El CHI EEH U2] for some specific cuts, where partial rates in semileptonic decays are 
expressed as weighted integrals over the B — > V s 7 photon spectrum, 



In the most recent analysis of this type, , the cut was chosen as P + < A for simplicity, and 
the weight function W(A, P + ) was calculated at complete two-loop order at the intermediate 
scale Hi. The second term on the right-hand side of this equation denotes a residual hadronic 
power correction (rhc), which was absorbed into the weight function in that reference. In the 
above relation T u denotes the partial semileptonic decay rate and (l/r s (E*)) (dT s /dP + ) the 
normalized photon spectrum in radiative decays, where P + = Mb — 2E 1) E 1 is the photon 
energy in the _B-meson rest frame, and the total decay rate T S (E*) is defined to include all 
events with > E* = rrib/20. It is beneficial to use the normalized photon spectrum instead 
of the absolute spectrum, because the weight function as defined in (JTJ) is independent of 
|K&Ksl) possesses a well-behaved perturbative expansion [T2|, and because the normalized 
photon spectrum can be determined with better accuracy than the absolute one [0]. In order 
to determine \V u b\ from relation (0) both T u and the normalized photon spectrum enter as 
experimental input, while the weight function W(A, P + ) and the residual hadronic correction 
r rnc are theoretical quantities. 

In this paper, we extend the technology developed in [T2] to derive the weight function for 
an arbitrary kinematic cut that includes events with P + ~ Aqcd and P_ 3> Aqcd- In general, 
the weight function W(A,P + ) (and the correction r r h c ), as well as the integration limit A 
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in (JTJ) depend on the particular cut that is used. For example, if we consider a cut on the 
charged-lepton energy Ei > E Q then A = Mb — 2E , or if a cut on hadronic invariant mass 
Mx < M Q is considered then one needs A = M . A quick way to obtain some intuition on 
how W(A, P + ) and A depend on the specific cut qualitatively is to consider the B — > X u l~v 
phase-space depicted in Figure [TJ Note that P = Mb in radiative decays, so that the phase- 
space of B — > Xg'-f consists only of the vertical boundary on the right-hand side of the plot. 
This picture explains correctly how the maximal P™ ax = A depends on the nature of the cut, 
and where to expect some kinks in W(A, P+), for example for a combined cut on Mx and q 2 . 

In the next section we will derive an expression for the weight function at leading power, 
which follows from exact factorization theorems for differential rates in inclusive B decays 
in the shape-function region. The formula for the weight function is thereby also valid to 
all orders in perturbation theory. We compute its explicit form to next-to-next-to leading 
order (NNLO) at the intermediate scale /i, and to next-to leading order (NLO) at the hard 
scale fih i n renormalization-group (RG) improved perturbation theory, including three-loop 
running effects. This approximation is sufficiently precise for phenomenological applications 
since effects at subleading power become as important as higher-order perturbative effects 
at leading power. In Section |3] we discuss first-order power corrections. Kinematical power 
corrections enter the weight function, while hadronic power corrections give rise the quantity 
r r hc in (HJ. We then apply our results to a few examples of kinematic cuts in Section |U 
and perform an analysis of theoretical uncertainties on \V u b\ by using a simple model for the 
experimental inputs. 

2 The weight function and factorization 

In this section we will adopt the leading-power approximation. The second term on the right- 
hand side of relation (JI} is then absent, because it collects contributions from subleading shape 
functions to both radiative and semileptonic decay rates, which start at order l/m b . 

2.1 Derivation 

We start the discussion by stating the exact factorization theorems for the fully differential 
leading-power decay rate in B — > X u l~v decays 
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and for the normalized B 



X s ^i photon spectrum 
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The superscript (0) indicates that these equations are valid to leading power. In (J5J) H u (y, e, hh) 
denotes the hard function, which collects all matching corrections at the scale Hh and depends 
on the kinematic variables 

P--P+ _ 2Ei 

V ~ M B -P + ' £ ~ 1 ~ M B -P + ' (4) 

for which the phase-space is < e < y < 1. The jet function J(p 2 ,fa) contains distributions 
that act on the shape function S(u, Hi) and depends on y via p 2 = ymf ) (P + — Cj) . In formula 
for the photon spectrum, on the other hand, the hard function is denoted by H r (E*, /i h ), and 
the argument of the jet function is p 2 = m b (P + — u), independent of y because P_ = M B in 
radiative decays. Three powers of the 6-quark mass appear due to phase-space integrations 
in the total rate T S (E*). Renormalization-group (RG) running effects between the scales Hh 
and Hi build up the functions U(Hh, Hi) and ar{Hh, Hi) hi the factorization formulas, which are 
such that U = 1 and a-p = in the limit fa ~ y t^h- Expressions for them will be given below. 

We start the derivation of the weight function by considering the factorized expression for 
an arbitrary partial differential decay rate 
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cut 



A Vm^[P+\ £ m ax[P+,J/] . . 

r r r/ 3 r ( } 

dP+ dy de " " , . (5) 

' ' dP + dyde v ; 







The integration limits y ma , x [P+] and s max [P + , y], as well as A depend on the specifics of the cut. 
(The use of squared brackets for these quantities is supposed to remind the reader that these 
functions are defined differently for different kinematic cuts.) We are now going to rewrite 
this semileptonic decay rate such that it resembles a weighted integral over the normalized 
photon spectrum. Clearly the biggest obstacle is that the argument of the jet function J in (J2J 
depends on the kinematic variable y, while it does not in (J3J). The solution to this problem has 
been presented in jT5j, where it was shown 1 that a factorization of the integrated jet function 
in the sense that 

/ rrib^l \ f n ( mdtt — k) \ 
J V /^~ + ^ V 1$ ' / 

can be achieved for arbitrary Q, if Y(k, Iny, //j) is allowed to be a distribution in the variable 
k. Here, 

f Q 2 \ f Q2 

j ln-=-,^i = / dp 2 J(p 2 ,fii) ■ (7) 



2 ' - 

K J Jo 

We will discuss the nature of Y(k, \ny, Hi) and its perturbative expansion to two-loop order 
in the next section. The strategy for the derivation of the weight function is to interchange 
integrations in so that the P + integration acts only on the jet function J, and we can 
make use of (jfij). Before this can be done, however, it is necessary to transform the other 

lr The quantity considered in that reference was f(k) — L dy dey~ 2ar H{y,e)Y{k 1 \ny). 
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P + -dependent terms in the partial decay rate into P + -independent ones. This can be achieved 
by inserting 



pA /-A pM B -P+ 

/ dP' 5{P' -P+)= / dP' / dq 5'(q + P'- M B) 

Jo Jo Jo 



into and by replacing P + — > P' in y max [P+], £max[P+] and in the kinematic prefactor 
(Mb — P + ) 5 . It is now an easy exercise to interchange the integrations 

A P+ M B -P+ A M B -Cb min(A,M s -g) 

dP+ I dou I dq...= I du I dq / dP+... . (9) 







Next, we apply (jOJ, interchange the k and P + integrations, and undo the steps in ©. Finally 
the integrations over q and P' can be carried out, which identifies P' = P + + k. As a result 
we arrive at an expression for the semileptonic decay rate, which is given as an integral over 
the product of the normalized photon spectrum in (jSJ) and 



A-P+ 





where 

!/ m ax[P'] e max [P',y] 

F(P , ,k, f JL h ,in) = (M B -P') s j dy J dey- 2a ^^H u (y } e,fi h )Y(kMy^ t ) . (11) 



The above two formulas enable us to calculate the weight function in an automated fashion. 
The procedure is as follows: first, the integration limits A, ?/ max [P + ], and £ max [P + ] are specified 
from the kinematics of the cut. For the next step it is helpful to decompose Y(k, hiy, /ij) = 
£\ T>i(k, Hi) Yi(\ny, //j), where T>i(k, Hi) are distributions in k and independent of y. (For 
example, we will see below that there are only three different distributions in the perturbative 
expansion of Y(k, \ny, fii) to two-loop order.) Likewise, we decompose F(P' , k, fih, Hi) = 
^2iT>i(k, /a) Fi(P', Hh, aO- It is then straight-forward to calculate the functions Fi(P',Hh,Hi) 
by integrating over e and y in (II ip . Finally the integration over k is performed in (|1U|). where 
the distributions Vi(k,Hi) act on Fi(P + + k,Hh,Hi)- 

Equation (jlOj) together with (jUJ) is the main result of this paper. We stress that both 
formulas are exact factorization theorems (in the convolution sense), valid to all orders in 
perturbation theory, and (J1(J|) gives the exact leading-power weight function for arbitrary 
cuts. 

2.2 Perturbative calculation 

The computation of the kernel Y(k,lny, h) in © requires the knowledge of the jet integral 
j(L,n)- The jet function J(p 2 ,Hi) has been calculated to one-loop order explicitly [HE], and 
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the non-constant part of j(L,fi) at two-loop order can be extracted from RG evolution jjj. 
Because the jet integral is of central importance to the present work, we find it legitimate to 
re-derive its dependence on L in detail. In particular, from the RG equations for the shape 
function and the leading-power current in soft-collinear effective theory (SCET, see also 
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(12) 



, 1 ^ + 2 7 V) + 2r c (/i)ln^ 



du m& J(m&(P+ — to), /x) S(uj, /x) , 



we can derive the RG equation governing j(L, /x) with L = In Q 2 / fi 2 and arbitrary Q 2 . Above, 
r c (/i) is the cusp anomalous dimension [Tlj, which has been calculated to three- loop order [To] , 
the remaining anomalous dimensions 7s(/x) and 7 / (/ i ) are known to two-loop order 
(For a definition of the star distribution that acts on the shape function see (llfij) below.) When 
combining the two RG equations, one finds 



d 



r 1 d z 

j(L^) = -2[T c (y)L + 1 J (y)]j(L,y)-2T c / — [j(L + ln(l - z),fi) - j(L, /x)] , (13) 

Jo z 



dln/x 

with 7 J = 7' — 75. This integro-differential equation can be solved perturbatively by choosing 
a polynomial ansatz for j(L,/x). As a result the integral in the above equation leads to the 
appearance of the Riemann zeta-function ( n . Specifically, after expanding the QCD /3-function 
and anomalous dimensions as 



da s (n) 
dlny 



n=0 V / n=0 ^ J 



n+1 



(14) 



and similarly 7 J (/x), we find to two-loop accuracy 



j(L,fj) = 1 + 



+ 



Air 



(15) 



^ ) +(^ 1) (7o J -/5o)+7 1 J 



7T" 

~6 



r 7o + C 3 rjj l 



+ \ [i J M - AO + ^r + r, - £rg) l 2 - - ^ J ) r ^ 3 + ^ 



Here, = Cf{7 — it 2 ) is the one-loop constant, and the two-loop constant b^' is currently 
unknown. To determine it a multi-loop calculation will be necessary. However, this constant 
does not enter the two-loop result for the kernel Y(k, In?/, /Xj). 

Next, we need to find an ansatz for Y(k, In y, /x) which satisfies (JBJ). At tree-level j(L, /x) = 1 
and therefore K(A;,lny,/x) = is also independent of y. It follows that the constant b^ 
cancels in the relation © to two-loop accuracy. Beyond the tree approximation the integrated 



■,(2) 
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Y(k, Iny, //) (defined equivalently to (J7|)) must pick up a logarithmic dependence on Q, as can 
be seen by interchanging the integrations in (JBJ); but Y(k, \ny, //) itself must not depend on Q. 
Objects that accomplish that are already known from the jet function J(p 2 , //), and are called 
"star-distributions" [TH] (see also jUGSE]). Their definitions are such that, when integrated 
over an interval Q, they act on a function <f)(k) as 



dk 



dk I — 



1 m fe /c 
t- m — — 
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0(0) 
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Therefore, the ansatz for the jet kernel reads (here with L = \ny for brevity) 



Y(k,L,fj) = 8{k) 



+ 



4:71 

4tt 



cS\L)6(k) + cY>(L)(i 



(!)/ 



c<i\L)6(k) + c?\L) 



(2), 



[M 2 /«ii,] 



[M 2 /mi,] 



(17) 



+ 4 2) (L) 



1 , m b k\ 
k 1*) 



[li 2 /m b ] 



In ffTUj) it acts on all /c-dependent objects, i.e., on the prefactor (M B — P + — k) 2 and on the 
integration limits y max [P + + k] and e max [P + + k,y]. To determine the coefficients c\ (L), 
however, we consult ©, where the star distributions act on the /c-dependent logarithms. We 
find 2 



c«(L) =r L, c£\L) = 
-T L(P - T L) 



.(2) 



7o L + y L > 



ci 2) (L)= [r 1 -/3 o7o J ]L + 
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J 7)( 1 ) R 



L + 



7T 
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L 2 + ^ 



7o J 



L 3 + —L A . 



At this point we interrupt the discussion briefly and consider the result for the weight function 
in (jlOJI again. As mentioned earlier, the 6-quark mass enters the normalization of the weight 
function since we also normalized the photon spectrum. In order to avoid the renormalon am- 
biguities of the pole scheme it is favorable to use a low-scale subtracted quark-mass definition, 
for which we adopt the shape-function scheme jH] . The two mass definitions are connected via 



m 



pole 



m b (/!*) + fjL : 



+ ... 



71 



(19) 



2 The expressions for the coefficients are generalizations of the corresponding findings in |12j . As mentioned 
in footnote ^ the connection is that the integrations over e and y were immediately performed, leading to the 
replacement of L n — > T n and subsequent devision by To in the notation of |12j . This was possible because for 
a pure cut on P + one has e max = V and y max = 1 independent of P+, unlike the general case considered here. 
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and we will refer to mp^/z*) at //* = 1.5 GeV as mj for brevity, for the rest of this paper. 
Apart from the normalization to the weight function, ra& also enters through radiative loga- 
rithms in the hard functions and through star distributions in the jet kernel. Since the hard 
functions are only kept to one-loop order, the above redefinition of m& does not affect the 
result. Y(k,lny, [Ai), on the other hand, is given to two-loop accuracy in (fT7J). It follows that 
cl 1 ^ receives a contribution at this level, 



C F a s (/j h ) fi* 
7o H J- o 



7T 



2 



(20) 



In the remainder of this section we collect the other ingredients of (fTU|) for completeness. 
The hard function in semileptonic decays is known to 1-loop order in perturbation theory and 
given by 



H u (y,e,n h ) = 12(y-e)(l-y + e)\l + 
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The hard function for the normalized photon spectrum reads 

C F a s (/i h ) 
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3 * 
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/is(^) 



(22) 



Here, 5* = 1 — 2£7*/m& = 0.9, and the functions fij(8*) capture effects from operator mixing 
[TU] and can be found in this notation in jB]. Note also the term proportional to a s (fih) n*/ m b, 
which ensures that the three powers of m;, in (jlOj) are defined in the shape-function scheme. 
Finally, the RG function a>r is defined as the integrated cusp anomalous dimension from In 
to ln/x/j, which yields 



In — ; — r H 
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(23) 
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When combining the various quantities into (JTUJ) the result should be re-expanded in a s to 
the order in which we are working. However, it is convenient to treat ar(/i/ l ,/ i i) as a running 
"physical" quantity (similar to a s (//)), which is not expanded. This is the same approach as 
put forward in [T2], and we will use it in the phenomenological applications in Section 0J 



Expansion coefficients of the anomalous dimensions. To three-loop order, the coeffi- 
cients of the (3 function read 
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C F + '— C A ) n) , 



where nj = 4 is the number of light flavors, Ca = 3 and Cp = 4/3. The cusp anomalous 
dimension to three- loop order is given by (14^ ITS] 



r = 4C F , 
r 2 = 16CV 
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We also need the anomalous dimension coefficients for the integrated jet function j(L, //), 
which are [S| 



7o = -3^f 
7i J = C, 



(26) 
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3 Power corrections 

Factorization theorems exist at each level of power counting for differential decay rates in inclu- 
sive heavy-quark decays. We differentiate between two types of power corrections, "kinemati- 
cal" and "hadronic" . The first class arises simply because we have restricted our discussion to 
a particular portion in phase-space, where P + <C Mb- These corrections are power suppressed 
in the shape-function region, but are of leading power when integrated over a domain that is 
comparable to Mb (OPE region). A different way of thinking about kinematical corrections 
is to view them as the equivalent of the factorization theorems (J2J) and Q with subleading 
hard and jet functions. However, no complete scale separation has been achieved for these 
power corrections yet, but the products of subleading hard and jet functions are known in 
fixed-order perturbation theory and to all powers from [T%1 119j. Since kinematical corrections 
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start at ot s (p) and are numerically small for all prominent cuts, this approximate treatment 
suffices. The scale p, is typically near a/^vV^cd, but independent of //j and 

The second class of corrections comes from subleading hadronic structure functions |2*T] 
122 1221 12H 123 EH 123 EH] ■ Already at first subleading order there are multiple independent 
shape functions entering the calculation of the differential decay rates. Furthermore, they 
appear in different linear combinations in the semileptonic and radiative cases, so that at 
this stage a weight function cannot relate semileptonic decay rates to radiative ones alone. 
In equation (JTJ) we have therefore added a second term labeled r rhc . Note that this term 
is different from the subleading shape-function contribution to the semileptonic decay rate 
(denoted r^ adr in jHj) since contributions from subleading hadronic corrections to the photon 
spectrum are convoluted with the leading-power weight function and must be subtracted. As 
a result the residual hadronic power corrections are collected in r rnc . Currently the hard and 
jet functions in the subleading shape-function contributions to the decay rates are only known 
at tree-level. 



3.1 Kinematical corrections 



Since kinematical corrections come with the leading shape function, the corresponding contri- 
bution to the weight function W km (A, P + ) can be calculated without the introduction of any 
hadronic uncertainty. These terms start at order a s (p,) and we only compute the correction 
to this order. It thus follows that W k[n <g> dvf ] = I^ in - W {0) ® dTf n , symbolically. (Here and 
below the abbreviation dT s in statements within the main text denotes the normalized photon 
spectrum.) The relevant expressions for the differential decay rates have been collected in [Hj 
and may be written as (here and below a r = ar(fai,fa) f° r brevity) 



cut 



N u I dP+ 



2/max[P+] £max[P+,2/] P+ 

dy J de y- 2ar (M B - P+) 4 J du K u (x, y, e) S{u) 



kin 



1 dT a 



N, 



(M B - P + 



mi 



du K s (x) S(lo] 



with x 



CO 



M B -P+ 



(27) 



where we have collected all terms in the functions K u (e,y,x) and K s (x), and abbreviated the 
different norms by N u and N s . The leading weight function at tree-level is taken from ([10)1 . 
We now transform the expression for r„ in — ® dT^ m in three steps: First, we interchange 
the order of the P + and Cj integrations. Next, the integration variable P + is substituted 
by k = P + — Cj. Finally the variable Cj is renamed by P + . After this has been done (no 
manipulation is needed for the term W km (g) dT^), we find 



W khl {A, P + 



mi 



A-P+ 



N s (M B - P + 



dk (M E 



kf 



dy 



de y 



-2a T 



X 



M R -P+ 



y,e)-Hi°\y,e)K s 



M f 



P,-k 



(28) 
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where Hu\y,e) denotes the tree- level part of the hard function in (|2T|). 

We are now going to restrict the calculation to the first-order power corrections because 
of mixing effects with hadronic power corrections at higher orders. For example, the second- 
order power correction to the weighted integral over the photon spectrum contains a term 
jykm(i) c /r^ adr ( 1 )_ (a s before, the superscript denotes the order in power counting.) For 
such a term we would require a compensation in the residual correction r r h c at order ot s (jx), 
which goes beyond the scope of this paper. Including only first-order power corrections is not 
a bad approximation; the studies in [Hj have shown that the full kinematical corrections can 
be approximated very well by including the first term in the power expansion only. At this 
level the functions K u (x,y,e) and K s (x) involve only two different functional dependences on 
x, a constant and a term proportional to In x. We find 



WKkin,l)( A) P 4 



G 2 F 



nr. 



C F a s {n) 



A-P+ 



192tt 3 (M B - P+) 3 4vr 



dk (M B -P+- k) 4 



o 



ymax I-* -f- tmax [P++k,y] 

x J dy J dey~ 2ar 
o o 



A(y,e) + B(y,e)ln 



(29) 



with 
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+ 4- [3 + 5e- (A + e)y + y 2 ] Iny 



+ 12(y-e)(l 
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48 



iCMCfAvh) 
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y 



i6jc^)p 



^(1 - y)(-3 + 5y - 5s) - |(1 - y + e)(y - e) [ ^% . 

y 3 [ c 7*(nh)\ 2 



In these expressions Cj(/i^) denote the (effective) Wilson coefficients of the relevant operators 
in the effective weak Hamiltonian for B — > X s 7 decay. Charm-loop penguin contributions to 
the hard function of the photon spectrum are encoded in the functions g±(z) and gi(z\ which 
depend on the variable z = (m c /mb) 2 ^Hj- They are 
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1 

,'/2( j / dx(l- x) 
'o 



x \zJ 2 



-2 arctan 2 y/t/(4 - t) 



G(t) 



2 In 



(v^+ v / t 3 4)/2 



in 



(30) 



t < 4 



-I ; t >4. 



This concludes the calculation of the weight function. 
3.2 Residual hadronic corrections 

There are four different hadronic structures entering the first-order power corrections to the 
differential decay rates at tree level. Following [21] we denote them by (A — u)S(uj), t(uj), 
u(uj), and v(u). The first one in this list involves the leading shape function and the heavy- 
quark parameter A. This parameter ensures that (A — uj)S(u) has a vanishing norm, as all 
subleading shape functions must have. It is possible to absorb the effect of this function into 
a weight-function contribution, which in turn only depends on (A — P+); however, because of 
the zero-norm constraint it is difficult 3 to assign the correct numerical value of A in that case. 
Instead, we keep this contribution together with the subleading shape functions in r r h c - A 
straight-forward calculation yields 



rhc 



GlUifih,^) f A 



cut 



16vr 3 



dP + {M B - P 4 



2/max[-P+] £max[P+:?/] 

) A j dy f dey-^y—± 
o o 



x { 2(1 - y)(y - e) (A - P + )S{P + ) + 2(1 - y) [ y - e + - ) t(P+) (31) 



e 



+ {l-y)(l-y + e)u{P + )+(l + y 



(2-y + y 2 ))v(P 



+ y(l-y + e)m 2 s S'(P + ) 



where we have kept the RG-evolution function y 2ar , as well as U(fih,f^i) [HI- Its structure is 
such that 

lnU(fi h ,fii) = 2S v (/i h , fa) - 2a v (/i h , ^) In 2ay(// fe ,^) , (32) 

and we evaluate to leading-order the Sudakov exponent (with r = a s (/ij)/a s (/z/i)) 



4p5 [a s {fi h ] V r 



^1 Pi \ /•, ,i \ , Pi i 2 

1 o Po / ^Po 



(33) 



3 In a specific default model for subleading shape functions was adopted such that this problem is 
avoided. A default model of this kind exists - but is different - for each kinematical cut. 
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and the RG function ay = r y' /(2/3 ) lnr analogous to ()23j). The use of the anomalous di- 
mension 7q = — 5CV of the leading SCET current is a good approximation because most 
of the subleading operators in SCET are build from the leading SCET current and thus 
share the same anomalous dimension. (A complete resummation requires knowledge of the 
anomalous dimension matrix of all subleading operators in SCET. For some discussion, see 
e.g. |23 EDI EH ■) 

We have included a small correction from a finite s-quark mass, leading to the expression 
proportional to m 2 s S\P + ). The appearance of hadronic structure functions in (j3Tj) introduces 
some irreducible uncertainties in phenomenological applications of our results. While in prin- 
ciple the leading shape function S(u) can be extracted from the photon spectrum, see e.g. jH], 
the forms of the shape functions t(u), u(uj), and v(uj) are unknown. What we know is that 
their norms vanish and their first moments are given in terms of the heavy-quark parame- 
ters n\ and A2. In order to estimate effects from higher moments we define functions h n (u), 
n G {t, u, v}, via 

i(u) = -A a S'{u>) + h t {u) , u{u) = -2& S'(u) + h u {u) , v(u) = A 2 S'{u) + h v (Qj) . (34) 

As long as each of the h n (Cj) have zero norm and first moment the constraints on the subleading 
shape functions are satisfied. 

Let us now adopt a default model for the terms in (|31|). For simplicity we use for S{uj) the 
familiar exponential-type functional form 

1 v / , &\ . , ,. 3A 2 

8 < u)= xw)[v exp r 6 iJ' mth (35) 

with parameters A and //|. By construction the moment constraints on (A — u>)S(u) are 
respected. The default model for (p^ljl is defined as h n (il)) = in all three cases. 

In order to estimate the uncertainty introduced by adopting a specific model we closely 
follow the procedure of H2], where four different functions hi(u), i = 1, ... ,4, were con- 
structed with vanishing norm and first moment. A variation of the functional form of the 
subleading shape functions was then achieved by setting h n (u) — > ±hi(u), where n G {t,u,v} 
and % G {1, . . . ,4}. Combinatorially this means we have 9 3 = 729 different models for the 
set of subleading shape functions. The estimator for the hadronic uncertainty is the maximal 
deviation from the default result when sampling over all models. 



4 Examples and Discussion 

Let us demonstrate the phenomenological implications of the result ([TO]). (J2SJ), and (|3*T|) by 
applying the main relation (0) for a few examples of typical kinematic cuts. To disentan- 
gle theoretical uncertainties from experimental ones we are going to pretend that both the 
normalized photon spectrum and the semileptonic partial decay rate were measured with no 
uncertainty. In particular we consider the photon spectrum given as 

1 dT s 1 b b J (P+\ b ~<~ 

r s (E*) dP + ~ A^r(6 7 ) \a^j 
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Figure 2: Examples of the weight function for different kinematic cuts. LEFT: Cutting on 
P + < A = 0.66 GeV and E\ > E . From top to bottom the four functions are for E = 0, 
E = 1 GeV, E = 2 GeV, and E = (M B - A)/2. RIGHT: Cutting on M x < M = 1.7 GeV, 
q 2 > <2q, and E\ > 1 GeV. The three functions are for q 2 . = (top), g 2 = 8 GeV 2 (middle), 
and ql = (M B - M ) 2 (bottom). 



with A 7 = 0.77 GeV and 6 7 = 2.5. This model describes the experimental data by BaBar [32| 
and Belle reasonably well. Furthermore we need as inputs the heavy-quark parameters 
A 2 = 0.12 GeV 2 , /4 = (0.25 ± 0.10) GeV 2 , and the quark masses m h = (4.61 ± 0.06) GeV 
[3IE1II1I, m s = (90 ± 25) MeV jSgEH], and m c /m h = 0.222 ± 0.027 0. Here m b and /i 2 are 
defined in the shape-function scheme at a scale = 1.5 GeV, while m s is evaluated in the 
MS scheme at 1.5 GeV. The ratio m c /mb is also evaluated in the MS scheme, where it is scale 
invariant. For the strong coupling a s (fi) we use three-loop running from a s (Mz) = 0.1187 
down to 4.25 GeV, apply matching corrections onto a 4-flavor theory, and then run to fih, f^i, or 
/2. Default values for these scales are taken to be mb/V2, 1.5 GeV, and 1.5 GeV, respectively 
To assign a perturbative error, the scales are varied around these default settings by factors 
of ^2 and 1/^/2. In all cases considered below, the analysis of uncertainties closely follows 
|12j . with similar outcome. 



Cutting on P+ and Ei. First, let us consider a cut on P + < A with A = 0.66 GeV together 
with a cut on E\ > E Q . In terms of the kinematic variables y and e this means that 

?/max[ J P+] = 1 , £max[i+, y] = Him ^1 - M ^ p ,VJ ■ (37) 

Using the central values of the input parameters and scales, the resulting weight functions are 
depicted on the left-hand side of Figure El for a few examples of E Q . There is an integrable 
singularity at the endpoint P + —>■ A ^TJ ^] if the lepton cut E is small. In the limit 
Eq — > (Mb — A)/2, corresponding to a pure cut on lepton energy, the weight function vanishes 
at the endpoint and the singularity disappears. 

We now study the case E = 1 GeV in more detail. This is a useful example because the 
partial branching fraction for this particular cut has already been measured [3*7] . First, we 
investigate the perturbative uncertainty of the right-hand side of (P), which is obtained by 
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studying the residual scale dependence entering via the weight function. (It is clear that the 
relevant quantity is the entire integral, and not the values of the weight function for individual 
values of P + .) We observe that the sum of the weighted integral over the photon spectrum 
W ® dT s = 44.58 ps -1 and the residual hadronic corrections r r h c = —5.65 ps -1 is very stable 
under scale variations, but not each of the two terms alone. The NNLO approximation for the 
weight function introduces roughly a 4% uncertainty. The error from the LO approximation 
to the power suppressed corrections is numerically of equal magnitude, but tends to cancel a 
large portion of the scale sensitivity. This leads us to interpret the perturbative error on the 
convolution integral as also the perturbative uncertainty on the sum of both terms, to avoid 
counting this error twice. The hadronic uncertainty on the residual term r r h c is obtained 
by taking the maximal deviation from the central value when sampling over a large set of 
models for the subleading shape functions, as outlined in Section l3~2l Finally we also vary the 
numerical value of and the remaining input parameters m s , m c , and fi 2 within their stated 
errors. This yields 

= 38.93 ^igg [pert.] ± 1.42 [hadr.] t\'Qj [m b ] 63 [P ars -1 ■ (38) 

P+ < 0.66 GeV, 
Ei > 1 GeV 

Further uncertainty enters in practice, because the photon spectrum cannot be measured over 
the entire range, but only over a certain window around the endpoint E 1 = Mb/2. The 
normalized photon spectrum is then obtained using theoretical information on the fraction of 
events that fall into this window. The current precision for these fractions is about 6% jS], 
which impacts W <S> dT s directly. If we further assume that the left-hand side of relation ([TJ 
was given with no experimental uncertainty, we can extract \V u i,\ with only a theoretical error. 
For example, let us take the central value Br(P + < 0.66 GeV, E x > 1 GeV) = 1.1 • 10" 3 [37], 
and dismiss the experimental error. Using the average lifetime tb = 1-60 ps of the B meson, 
and taking the normalization uncertainty on the photon spectrum into account, we find 

\V ub \ = (4.20 ^ [theory]) • 10- 3 . (39) 



w ® dv s + r rhc 

ps -1 



Cutting on M x , q 2 , and E t . When cutting on M x < M , q 2 > q 2 , and E\ > E , the 
maximal value of P + is given by A = min(Mo, Mb — %/*?o > ~ 2E ). The phase space is 
such that 

V^[P + \ = min (l - Jp^tf, _ + p+) ), (40) 

ew*x[P+, y] = min ^1 - M ^ p » V 

The right-hand side of Figure 121 shows three different weight functions. In all cases M = 1-7 
GeV, close to the optimal value Mp, and E — 1 GeV as before. From top to bottom the three 
curves are for q 2 = (pure Mx cut), q 2 = 8 GeV 2 (mixed cut), and q 2 = (Mb — M ) 2 (pure q 2 
cut). Note that the integrable singularity of the pure P + cut in the previous discussion is gone, 
because this point is no longer the endpoint. Instead, the weight function has a kink, which 
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is expected from considerations of the phase space depicted in Figure [T] The endpoint A is 
much larger in this case, and the weight function vanishes there, so that we do not encounter 
any singularity anymore. 

As a second example of a \ V u b\ determination we consider the case M = 1.7 GeV, q 2 = 0, 
£o = l GeV, i.e., the top curve in the plot on the right-hand side of Figure 121 Because of 
phase-space restrictions the weight function has a kink at P + « 0.55 GeV. For the central 
values and default models we find W ® dT s + T rhc = (49.77 — 1.63) ps -1 , and the analysis of 
uncertainties yields 



W ® dT s + r 



rhc 



ps -1 



= 48.14 till [Pert.] ± 0.47 [hadr.] t\f 8 H tlf 9 [pars.] . (41) 

M x < 1.7 GeV, 
E l > 1 GeV 



Again, we must also add a 6% uncertainty to the norm of the photon spectrum. From the 
input Br[M x < 1.7 GeV, E x > 1 GeV) = 1.24 ■ 10~ 3 jSj follows 

\V ub \ = (4.01 t° il [theory]) • 10~ 3 . (42) 

The third and last example is the combined cut M = 1.7 GeV, q 2 = 8 GeV 2 , E Q = 1 GeV, 
whose weight function is shown as the curve in the middle of the right plot in Figure El In 
analogy to the previous cases we obtain W <8> dT s + T r ^ c = (28.18 — 4.97) ps -1 and 



w ® dr s + r 



rhc 



ps 1 



M x < 1.7 GeV, 
q 2 > 8 GeV 2 , 
E t > 1 GeV 



23.21 [pert.] ± 0.43 [hadr.] ^{ ® [m b ] t-Q^Q [P ars -1 ■ (43) 



Taking Br(M x < 1.7 GeV, q 2 > 8 GeV 2 , E t > 1 GeV) = 8.56 • 10~ 4 , which is the mean of the 
central values of the measurements reported in (S3 EH], leads to a larger value, namely 

\V ub \ = (4.79 [theory]) ■ 10- 3 . (44) 

It is of course understood that these values change once an analysis with the full exper- 
imental data is conducted. We believe, however, that a theoretical error of about ±5% is 
realistic and consistent with the previous work in [SI ^] . 



Concluding remarks. In the last section we have given results for the weight function 
W(A, P + ) for several commonly employed kinematic cuts. In all cases the charged-lepton 
energy was bound to exceed 1 GeV, which is typically used to identify semileptonic B decays 
in practice. We stress that for an additional cut on P + the photon spectrum is required over 
only a small window, where already precise data exists. For a cut on the hadronic invariant 
mass, on the other hand, the photon spectrum is also needed in a regime where its measurement 
is very difficult. Cutting away further events in the low P + region by virtue of an additional 
restriction on the leptonic invariant mass q 2 worsens the situation even more, because the 
relative importance of the high P + region is enhanced. It is therefore expected that the first 
example, the cut on P + , will ultimately lead to the most precise determination of IV^. 



16 



In summary, we have presented a formula in (|1U|) which is based on exact factorization 
theorems for the differential decay rates, and that allows for the calculation of weight functions 
for arbitrary kinematic cuts. Quantities entering this formula were evaluated with one-loop 
precision at the hard scale, and complete two-loop precision at the intermediate scale, including 
three-loop running effects in renormalizat ion-group improved perturbation theory. To achieve 
further precision we have also included first-order power corrections resulting from subleading 
kinematical and hadronic contributions. The details of the cut are encoded in three kinematic 
quantities, £ max [P+,y], ymax[P+], an d A. Once they are specified, one only needs to carry out 
a few integrations that lead directly to the weight function. 

The use of relations such as ((TJ) circumvents the necessity for fitting models of the leading 
shape function to the B — > X s ^ photon spectrum and allows for the determination of |T4&| in 
a model-independent way at leading power. The quest for a precision measurement of \V u b\ 
requires a variety of different approaches. The results of this paper represent an alternative 
route to direct theoretical predictions of partial decay rates. 
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